N = [2:65,70:5:200];

s = zeros(1,length(N));

for i = 1:length(N)
    fileName = strcat('traps',num2str(N(i)),'.dat');
    [theta, phi] = loadTrapData(fileName);
    [x,y,z] = fromSphToCart(1, theta, phi);
    
    for j = 1:length(x)
        for k = j+1:length(x)
            s(i) = s(i) + log(2+eucDistanceCart(x(j),y(j),z(j),x(k),y(k),z(k)));
        end
    end
end

plot(N,s)

